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Abstract. By applying an idea of Borodin and Olshanski (2007), we study various 
scaling limits of determinantal point processes with trace class projection kernels 
given by spectral projections of selfadjoint Sturm-Liouville operators. Instead of 
studying the convergence of the kernels as functions, the method directly addresses 
the strong convergence of the induced integral operators. We show that, for this 
notion of convergence, the Dyson, Airy, and Bessel kernels are universal in the 
bulk, soft-edge, and hard-edge scaling limits. This result allows us to give a short 
and unified derivation of the known formulae for the scaling limits of the classical 
unitary random matrix ensembles (GUE, LUE/Wishart, JUE/MANOVA). 



1. Introduction 

We consider determinantal point processes on an interval A — (a, b) with trace 
class projection kernel 

n-1 

K n (x,y) = £ <fy(x)<fy(y), (1.1) 

j=0 

where cpo, <p\, ■ ■ ■ , <p n -i are orthonormal in L 2 (A); each <pj may have some depen- 
dence on n that we suppress from the notation. We recall (see, e.g., Anderson, 
Guionnet and Zeitouni 2010, §4.2) that for such processes the joint probability 
density of the n points is given by 

1 " 

p n (xi,...,x„) = — det K n (xi,Xj), 
n\ i,j=i 

the mean counting probability is given by the density 

p n (x) = n~ l K„(x,x), 

and the gap probabilities are given, by the inclusion-exclusion principle, in terms 
of a Fredholm determinant, namely 

£„(/) = F({ Xl x n } n / = 0) = det(I - IjKnlj). 

The various scaling limits are usually derived from a pointwise convergence of the 
kernel function K n (x,y) obtained from considering the large n asymptotic of the 
eigenf unctions <pj, which can be technically very involved. 1 Borodin and Olshan- 
ski (2007) suggested, for discrete point processes, a different, conceptually and 
technically much simpler approach based on selfadjoint difference operators. We 

2010 Mathematics Subject Classification. 15B52, 34B24, 33C45. 

x Based on the two-scale Plancherel-Rotach asymptotic of classical orthogonal polynomials or, 
methodically more general, on the asymptotic of Riemann-Hilbert problems; see, e.g., Deift (1999), 
Johnstone (2001), Collins (2005), Johnstone (2008), Anderson et al. (2010), and Kuijlaars (2011). 
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will show that their method, generalized to selfadjoint Sturm-Liouville operators, 
allows us to give a short and unified derivation of the various scaling limits for 
the unitary random matrix ensembles (GUE, LUE/Wishart, JUE/MANOVA) that 
are based on the classical orthogonal polynomials (Hermite, Laguerre, Jacobi). 

The Borodin-Olshanski Method. The method proceeds along three steps: First, 
we identify the induced integral operator K n as the spectral projection 

K n = l^^o^Ln) 

of some selfadjoint ordinary differential operator L„ on L 2 (A) . Any scaling of the 
point process by x = a n ^ + yi n (c„ ^ 0) yields, in turn, the scaled objects 

£„(/) = det(J - K n = 1 { _ 00<0) (L„), 

where L n is a selfadjoint differential operator on L 2 (A n ), A n = (a n ,b n ). 

Second, if A„ C A = (a, b) with a n — > a, b n — > b, we aim for a selfadjoint 
operator L on L 2 (A) with a core C such that eventually C C D(L n ) and 

L n u — > Lu (u s C). (1.2) 

The point is that, if the test functions from C are particularly nice, such a con- 
vergence is just a simple consequence of the locally uniform convergence of the 
coefficients of the differential operators L„ — a convergence that is, typically, an easy 
calculus exercise. Now, given (1.2), the concept of strong resolvent convergence (see 
Theorem 2) immediately yields, 2 if ^ Cpp(L), 

K "^A„ = 1 (-oo,0)(^) 1 A„ > 1 (-oo,0)(^)- 

Third, we take an interval / C A, eventually satisfying / C A„, such that the 
operator q) (£)!/ is trace class with kernel K(x, y) (which can be obtained from 
the generalized eigenfunction expansion of L, see §A.2). Then, we immediately 
get the strong convergence 

K„lj Klj. 

Remark 1.1. Tao (2012, §3.3) sketches the Borodin-Olshanski method, applied to 
the bulk and edge scaling of GUE, as a heuristic device. Because of the microlocal 
methods that he uses to calculate the projection l(_oo,o) (£)/ he puts his sketch under 
the headline "The Dyson and Airy kernels of GUE via semiclassical analysis". 

Scaling Limits and Other Modes of Convergence. Given that one just has to 
establish the convergence of the coefficients of a differential operator (instead of an 
asymptotic of its eigenfunctions), the Borodin-Olshanski method is an extremely 
simple device to determine all the scalings x — o~ n £ + }i n that would yield some 
meaningful limit K n lj — > Klj, namely in the strong operator topology. Other 
modes of convergence have been studied in the literature, ranging from some 
weak convergence of A: -point correlation functions over pointwise convergence of 
the kernel functions to the convergence of gap probabilities, that is, 

E n (l) = det(I - IjKnlj) -> det(7 - ljKlj) = £(/). 

From a probabilistic point of view, the latter convergence is of particular interest 
and has been shown in at least three ways: 



2 By we denote the strong convergence of operators acting on L 2 . 
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(1) By Hadamard's inequality, convergence of the determinants follows di- 
rectly from the locally uniform convergence of the kernels K n (Anderson 
et al. 2010, Lemma 34.5) and, for unbounded /, from further large devia- 
tion estimates (Anderson et al. 2010, Lemma 3.3.2). This way, the limit gap 
probabilities in the bulk and soft edge scaling limit of GUE can rigorously 
be established (see, e.g., Anderson et al. 2010, §§3.5 and 3.7). 

(2) Since A 1— >■ det(J — A) is continuous with respect to the trace class norm 
(Simon 2005, Thm. 3.4), K n lj — > Klj in trace class norm would generally 
suffice. Such a convergence can be proved by factorizing the trace class op- 
erators into Hilbert-Schmidt operators and obtaining the L 2 -convergence 
of the factorized kernels once more from locally uniform convergence, see 
the work of Johnstone on the scaling limits of the LUE/Wishart ensembles 
(2001) and that on the limits of the JUE/MANOVA ensembles (2008). 

(3) Since ljK n tj and IjKIj are selfadjoint and positive semi-definite, yet 
another way is by observing that the convergence K n lj — »■ Klj in trace 
class norm is, for continuous kernels, equivalent (Simon 2005, Thm. 2.20) 
to the combination of both, the convergence K n tj — > Klj in the weak 
operator topology and the convergence of the traces 



Once again, these convergences follow from locally uniform convergence 
of the kernels; see Deift (1999, §8.1) for an application of this method to 
the bulk scaling limit of GUE. 

Since convergence in the strong operator topology implies convergence in the 
weak one, the Borodin-Olshanski method would thus establish the convergence of 
gap probabilities if we were only able to show condition (T) by some additional, 
similarly short and simply argument. Note that, by the ideal property of the trace 
class, condition (T) implies the same condition for all J' C /. We fall, however, 
short of conceiving a proof strategy for condition (T) that would be independent of 
all the laborious proofs of locally uniform convergence of the kernels. 

Remark 1.2. Contrary to the discrete case originally considered by Borodin and 
Olshanski, it is also not immediate to infer from the strong convergence of the 
induced integral operators the pointwise convergence of the kernels, 



to prove a limit law p n {i) dt — > p(t) dt for the mean counting probability. Using 
mollified Dirac deltas, pointwise convergence would generally follow, for continu- 
ously differentiable K n (^,rj), if were able to bound, locally uniform, the gradient of 
K n (^,n). Then, by dominated convergence, criterion (T) would already be satisfied 
if we established an integrable bound of K n (£, on /. Since the scalings laws 
are, however, maneuvering just at the edge between trivial cases (i.e., zero limits) 
and divergent cases, it is conceivable that a proof of such bounds might not be 
significantly simpler than a proof of convergence of the gap probabilities itself. 




(T) 



For instance, in §2, we will need just a single such instance, namely 

X„(0,0) ->K(0,0), 



(Ko) 
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The Main Result. Using the Borodin-Olshanski method, we will prove the fol- 
lowing general result on self adjoint Sturm-Liouville operators; a result that adds 
a further class of problems to the universality (Kuijlaars 2011) of the Dyson, Airy, 
and Bessel kernel in the bulk, soft-edge, and hard-edge scaling limits. 

Theorem 1. Consider a selfadjoint realization L n of the formally selfadjoint Sturm- 
Liouville operator 3 

d 



Ln = -TxV {x) Tx) +qn[x) - Xn 

on A = (-00,00), A = (0,oo), or A — (0,1) with p,q n G C°°(A) and p(x) > for all 
x 6 A such that, for t 6 A and n — > 00, 

n- 2K '\ n ~u), n- 2K 'q n (n K t) ~ q{t), n 2K "p{n K t) ~ p(f) > 0, (1.3) 

asymptotically up to an error 0(n _1 ); f/ie scaling exponents should satisfy 

k + k' + k" = \. (1.4) 

We assume that these expansions can be differentiated 4 at least twice, that the roots of 
q(t) — co are simple and that the spectral projection K n = l^^g) (L n ) satisfies trK n = n. 
Let x = a n ^ + y„ induce the transformed projection K n and let x = n K t induce the 
transformed mean counting probability density p n . Then the following scaling limits hold. 

• Bulk scaling limit: given t 6 A with q(t) < cv, the scaling 

n K-l I 

]i n = n K t, cr n = -^jj, p(t) = 

yields, for any bounded interval J, the strong operator limit K n lj — % i^Dysonl/- 
Under condition (Ko) and if p has unit mass on A, there is the limit law 

p„{t)dt p{t)dt. 

• Soft-edge scaling limit: given f* £ A with q(U) = co, the scaling 




_ 2 

u n x = n K t*, a„ = n K ~3 



P(t*) 



1/3 



yields, for s > —00 and any interval } C (s, 00), the limit K n lj — > K^ry 1/. 
Hard-edge scaling limit: given that A = (0, 00) or A = (0, 1) with 

p(0)=0, p'(0)>0, q n (x) =q(x) =y 2 x- 1 +0(1) (x ^0), (1.5) 
ffae scaling 

n P'(°) 2 T 



4am 2 *" VpW 
yields, 5 for any bounded interval fc (0,<x>), the limit K„lj — ^ Vessel 

^Srnce, in this paper, we always consider a particular selfadjoint realization of a formal differential 
operator, we will use the same letter to denote both. 

4 We say that an expansion f n (t) -/(f) = 0(1/ n) can be differentiated if /£(f) -/'(f) = 0(l/u). 
^Here, if ^ a < 1, the selfadjoint realization L„ is defined by means of the boundary condition 

2xh'(x) - «k(x) = o(jr" /2 ) (x-+0). (1.6) 



ON THE SCALING LIMITS OF DETERMINANTAL POINT PROCESSES 



5 



Remark 1.3. Whether the intervals / can be chosen unbounded or not depends 
on whether the limit operator Klj would then be trace class or not, see the 
explicit formulae for the traces in the appendix: only in the former case we get the 
representation of the scaling limit in terms of a particular integral kernel. Note 
that we can never use / = A since tr K n = n — > 00. 

Outline of the paper. The proof of Theorem 1 is subject of §2. In §3 we apply it to 
the classical orthogonal polynomials, which yields a short and unified derivation of 
the known formulae for the scaling limits for the classical unitary random matrix 
ensembles (GUE, LUE/Wishart, JUE/MANOVA). In fact, by a result of Tricomi, 
the only input needed is the weight function w of the orthogonal polynomials; 
from there one gets in a purely formula based fashion (by simple manipulations 
which can easily be coded in any computer algebra system), first, to the coefficients 
p and q n as well as to the eigenvalues k n of the Sturm-Liouville operator L n and 
next, by applying Theorem 1, to the particular scalings limits. 

To emphasize that our main result and its application is largely independent 
of concretely identifying the kernel K of the spectral projection l(_ 00/ o) (!<)/ we 
postpone this identification to the appendix: there, using generalized eigenfunction 
expansions, we calculate the IDyson, Airy, and Bessel kernels directly from the 
limit differential operators L. 

2. Proof of the Main Result for Sturm-Liouville Operators 

To begin with, we first state the expansions of the coefficients, which also 
motivate the particular regularity assumptions made in Theorem 1. 

Expansions. Since L n is a selfadjoint realization of 



with p,q„ 6 C°°(A) and p(x) > for x G A, we have C~(A) C D(L„). The scaling 



(with (6Aa parameter) yields A = (—00,00) in each of the following cases: 

• A = (—00,00); 

• A = (0,oo), k ^ 0, cr n -> 0; 

• A= (0,1), k = 0, a n -> 0. 

By virtue of t]) drj = K n (x,y) dy (with y — cr n rj + n K t), such a scaling induces 




x = cr„£ + n K t 



(2.1) 



a n K n (cr n ^ + n K t,cr n ri + n K t). 



Because of 



> and p{n K t) > 0, we see that K n = l(_ 00/ o)(£n) with 
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Now, if cr n = o(n K ~ 1 ^ 2 ), we get from (1.3) and (1.4) by Taylor expansion that 6 

L„u(£) = - u "(^) + a |^( ? (f)-a; + c7„n-Y(f)^)w(0 + o(l)- 
Bulk Scaling Limit. If q(t) 7^ to, by choosing 

<r n = cr n (t) = Tin*- 1 J 

we get L„u(f) = — «"(£) — S7r 2 «(£) + o(l); that is, L n w — >■ Lm with the limit 

d 2 , 

L = - ^ ~ S7r ' s = sign(o; - <?(£)). 

We recall that L is essentially selfadjoint on Cq°(A) and that its unique selfadjoint 
extension has absolutely continuous spectrum: cr(L) = u ac (L) = [— szr 2 , 00). Thus, 
for s = —1, the spectral projection %r OOi0 \(L) is zero. For s = 1, the spectral 
projection can be calculated by a generalized eigenfunction expansion, yielding 
the Dyson kernel (A. 3); see Lemma A.i. 

Limit Law. The result for the bulk scaling limit allows, in passing, to calculate a 
limit law of the mean counting probability p n (x) dx. We observe that, by virtue of 
p n {x) dx = p n {t) dt, the scaling x = n K t transforms p n (x) into 7 



p n (t) = n K p n (n K t) = n*- x K n (nH,nH) = ^£„(0,0) = -J ^ f (f }l £„(0,0). 

a n (t) Tz]j p(t) 

Thus, to get to a limit, we have to assume condition (Ko), that is, 8 
X„(0,0) ->• [q(t) < to] K Dyson (0,0) = [q(t) < to]. 
We then get, for q(t) 7^ to (that is, almost everywhere in A), 



Hence, by Helly's selection theorem, the probability measure p n {t) dt converges 
vaguely to pit) dt, which is, in general, just a sub-probability measure. If, however, 
it is checked that p(t) dt has unit mass, the convergence is weak. 

Soft-Edge Scaling Limit. If q{t*) = to, by choosing 9 

1/3 



we get L„u(£) = — u"(£) + + o(1); that is, L„m — > Lu with the limit 

^ 2 



r2 1 " f ■ 



6 If cr n = o(h k_1 ' /2 ), the error term is 0(max(cr„tt~' c / t7" 2 j; 1 ~ 2 ' f )), which amounts for 0(n _1 ) in the 
bulk scaling limit and for 0(n -1 ' 3 ) in the soft-edge scaling limit. 

7lt is the cancellation of the powers of n in the last equality, which called for assumption (1.4). 
^Following Knuth, the bracketed notation [S] stands for 1 if the statement S is true, otherwise. 
9 Note that, by the assumption made on the simplicity of the roots of 17(f) — w, we have ^'(f*) 7^ 0. 
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We recall that L is essentially selfadjoint on C^°(A) and that its unique selfadjoint 
extension has absolutely continuous spectrum: c(L) = c r ac (L) = (—00,00). The 
spectral projection l(_ 00/ o)(^') can straightforwardly be calculated by a generalized 
eigenfunction expansion, yielding the Airy kernel (A.4); see Lemma A.2. 

Hard-Edge Scaling Limit. In the case a = 0, we take the scaling 

x = cr n £ (y„ > 0), 

with c n =o(l) appropriately chosen, to explore the vicinity of this "hard edge"; 
note that such a scaling yields A = (0, 00). To simplify we make the assumptions 
stated in (1.5). We see that K n — l(_oo,0)(£n) with 

r 4 d ( . d \ Acr„ , , , . 

u = -jW- n d-i V^rd + W) {q{an0 ~ Xnh 

If we choose 

a ' 1 ~ n 4a; ' 

we get 10 L n u(£) = -4(^«'(gr))' + (^r 1 " 1) »(£) +o(l); that is, L„w ->■ Lw, 

We recall that, if a ^ 1, the limit L is essentially selfadjoint on C~(A) and that 
the spectrum of its unique selfadjoint extension is absolutely continuous: c(L) = 
C ac (X) = [— l,oo). The spectral projection can be calculated by a generalized 
eigenfunction expansion, yielding the Bessel kernel (A. 5); see Lemma A. 3. 

Remark 2.1. The theorem also holds in the case ^ a < 1 if the particular selfadjoint 
realization L„ is defined by the boundary condition (1.6), see Remark A.i. 

3. Projection Kernels Associated to Classical Orthogonal Polynomials 

In this section we apply Theorem 1 to the kernels associated with the classical 
orthogonal polynomials, that is, the Hermite, Laguerre, and Jacobi polynomials. 
In random matrix theory, the thus induced determinantal processes are modeled 
by the spectra of the Gaussian Unitary Ensemble (GUE), the Laguerre Unitary 
Ensemble (LUE) or Wishart ensemble, and the Jacobi Unitary Ensemble QUE) or 
MANOVA 11 ensemble. 

To prepare the study of the individual cases, we first discuss their common 
structure. Let P n (x) be the sequence of classical orthogonal polynomials belonging 
to the weight function zv(x) on the interval (a, b). We normalize P„ (x) such that 
(<pn><pn) — 1/ where cp n (x) = w(x) l ^ 2 P n (x). The functions <p n form a complete 
orthogonal set in L 2 (a, b); conceptual proofs of the completeness can be found, e.g., 
in Andrews, Askey and Roy (1999) (§5-7 for the Jacobi polynomials, §6.5 for the 
Hermite and Laguerre polynomials). By a result of Tricomi (see Erdelyi, Magnus, 
Oberhettinger and Tricomi 1953, §10.7), the P n (x) satisfy the eigenvalue problem 

~w(xjlx (pW^W^^W) = A nPn{x), X„ = -n(/ + \{n + l)f), 



I0 The error term is 0(n mm ( 1 > 2K )), which amounts for 0(n 1 ) if k' ^ 1 (as throughout §3). 



MANOVA = multivariate analysis of variance 
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where p(x) is a quadratic polynomial 12 and r(x) a linear polynomial such that 

w'(x) r(x) 
w(x) p{x)' 

In terms of cp n , a simple calculation shows that 

d ( d \ r(x) 2 r'(x) 

~dx V^dx^W) +^( X )^"( X ) = X n<pn{x), q(x) = + 

Therefore, by the completeness of the (p n , the formally self adjoint Sturm-Liouville 
operator L = ~f^p{ x )j^ + <?(*) has a particular selfadjoint realization (which we 
continue to denote by the letter L) with spectrum 

cr(L) = {A 0/ A 1/ A 2 ,...} 

and corresponding eigenfunctions cp n . Hence, the projection kernel (1.1) induces 
an integral operator K n with tr K n = n that satisfies 

K n = l(_oo,0)(^'n); L„ = L — A n - 

Note that this relation remains true if we choose to make some parameters of the 
weight w (and, therefore, of the functions <pj) to depend on n. For the scaling limits 
of K n , we are now in the realm of Theorem 1: given the weight zv(x) as the only 
input, all the other quantities can now simply be obtained by routine calculations. 

Hermite Polynomials. The weight is w(x) = e~ x on A = (—00,00); hence 

p(x) = 1, r(x) = —2x, q(x) = x 2 — 1, A„ = 2n, 
and, therefore, 

k = k' = \, k" = 0, p(t) = 1, q{t) = t 2 , co = 2. 

Theorem 1 is applicable and we directly read off the following well-known scaling 
limits of the GUE (see, e.g., Anderson et al. 2010, Chap. 3): 

• bulk scaling limit: if — \[2 < t < \/2, the transformation 

x = n l +n l ' 2 t 

induces K n with a strong limit given by the Dyson kernel; 

• limit law: the transformation x = n x l 2 t induces p n with a weak limit given 
by the Wigner semicircle law 

• soft-edge scaling limit: the transformation 

x = ±(2- 1/2 n~ 1/6 £+V2n~) 
induces K n with a strong limit given by the Airy kernel. 



2 With the sign chosen that p(x) > for x £ («,&). 
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Laguerre Polynomials. The weight is w(x) = x a e~ x on A = (0,oo); hence 

p(x) = x, r(x) =a.-x, q(x) = — — ^ ~, A„ = n. 

In random matrix theory, the corresponding determinantal point process is mod- 
eled by the spectra of complex n x n Wishart matrices with a dimension parameter 
m ^ n; the Laguerre parameter a is then given by a. = m — n ^ 0. Of particular 
interest in statistics (Johnstone 2001) is the simultaneous limit »!,n->co with 

m 

— >e>i, 

n 

for which we get 

/I // 1 / \ , s (0 - 1 - f) 2 

K = l, k> = -, k" = --, p(t) = t, q(t) = - — 4^-^-/ w = l. 

Note that 

^"^"H f ± = (v^±i) 2 . 

Theorem 1 is applicable and we directly read off the following well-known scaling 
limits of the Wishart ensemble (Johnstone 2001): 



• bulk scaling limit: if f_ < t < f +, 

27rt£ 



V(t+-t)(*-t- 

induces K n with a strong limit given by the Dyson kernel; 

limit law: the scaling x = nt induces p n with a weak limit given by the 
Marcenko-Pastur law 

1 



• soft-edge scaling limit: with signs chosen consistently as either + or — , 
x = ±n 1/3 0- 1/6 £ /3 £ + nt± (3.1) 
induces X„ with a strong limit given by the Airy kernel. 
Remark 3.1. The scaling (3.1) is better known in the asymptotically equivalent form 



x = a^ + }i, u = (y/m± V n ) 2 , °~ = W m ^ v 7 ") ( —7= ± 



1 1 



1/3 



which is obtained from (3.1) by replacing 6 with m/n (see Johnstone 2001, p. 305). 

In the case = 1, which implies f_ = 0, the lower soft-edge scaling (3.1) breaks 
down and has to be replaced by a scaling at the hard edge: 

• hard-edge scaling limit: if a. = m — n is a constant/ 3 x = £/(4tt) induces 

K n with a strong limit given by the Bessel kernel Cesser 

I3 By Remark 2.1, there is no need to restrict ourselves to a ^ 1: since </>>i(-y) = i*j;«(a;) with XhM 
extending smoothly to x = 0, we have, for a ^ 0, 

x" /2 (2x<p' n (x) - Kfi,(x)) = 2x 1+ei x' n (x) = 0(x) (x -> 0). 

Hence, the selfadjornt realization L n is compatible with the boundary condition (1.6). 
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Jacobi Polynomials. The weight is w{x) = x a (l — x)P on A = (0, 1); hence 

f \ /- \ / \ , „, / \ ice - (a + R)x) 2 a + B 
p(x) = x(l-x), r{x)=OL-{oc + p)x, q{x) = — _ ^ — , 

and 

X„ = n{n + ct + p + 

In random matrix theory, the corresponding determinantal point process is mod- 
eled by the spectra of complex n x n MANOVA matrices with dimension param- 
eters mi, mi ^ n; the Jacobi parameters a, f> are then given by a = m\ — n ^ 
and j6 = mi — n ^ 0. Of particular interest in statistics (Johnstone 2008) is the 
simultaneous limit m.\,mi,n — > 00 with 

— ^ >ee(o,i), —J > re (0,1/2], 

m\ + mi m\ + mi 

for which we get 

n n 1 + ,s * — / \ (6-T- (1 -2t)£) 2 1-t 

Note that 



Theorem 1 is applicable and we directly read off the following (less well-known) 
scaling limits of the MANOVA ensemble (Collins 2005, Johnstone 2008): 

• bulk scaling limit: if f_ < t < f +, 

27TTf (1 - t) £ , , 



nV(*+ -*)(*-*-) 
induces K n with a strong limit given by the Dyson kernel; 

limit law: (because of k = there is no transformation here) p n has a weak 
limit given by the law (Wachter 1980) 





soft-edge scaling limit: with signs chosen consistently as either + or 

-2/3 Crt±(l-f±)) 2 ' 3 

(t0(1-t)(1-0)) 
induces X„ with a strong limit given by the Azry kernel. 



Remark 3.2. Johnstone (2008, p. 2651) gives the soft-edge scaling in terms of a 
trigonometric parametrization of 6 and t. By putting 

= sm z |, t = sin" |-, 

we immediately get 

f± = sin 2 
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and (3.2) becomes 

x = ±a ± Z + t ±/ <r± = n~ 2/3 , 7, p . 

\ 4 sin q> sin if> / 

In the case 9 = T = 1/2, which is equivalent to m\ln,miln — > 1, we have 
£_ = and f+ = 1. Hence, the lower and the upper soft-edge scaling (3.2) break 
down and have to be replaced by a scaling at the hard edges: 

• hard-edge scaling limit: if a = m\ — n, /3 = mi — n are constants/ 4 x = 

£/ (4n 2 ) induces K n with a strong limit given by the Bessel kernel K^ sse {, by 

symmetry, the Bessel kernel K^se! is obtained for x = 1 — £/ (4n 2 ). 



A. Appendices 

A.i. Generalized Strong Convergence. The notion of strong resolvent convergence 
(Weidmann 1980, §9.3) links the convergence of differential operators, pointwise 
for an appropriate class of smooth test functions, to the strong convergence of 
their spectral projections. We recall a slight generalization of that concept, which 
allows the underlying Hilbert space to vary. 

Specifically we consider, on an interval (a, b) (not necessarily bounded) and on 
a sequence of subintervals (a n ,b n ) c (a, b) with a n — > a and b n — » b, selfadjoint 
operators 

L : D(L) C L 2 {a,b) L 2 {a,b), L n : D(L„) C L 2 (a n ,b„) -> L 2 (a„,&„). 

By means of the natural embedding (that is, extension by zero) we take L 2 {a n ,b n ) C 
L 2 (a,b); the multiplication operator induced by the characteristic function h a „,b„), 
which we continue to denoted by lj fln ;,„), constitutes the orthogonal projection 
of L 2 (fl, b) onto L 2 (a„, b n ). Following Stolz and Weidmann (1993, §2), we say that 
L n converges to L in the sense of generalized strong convergence (gsc), if for some 
z 6 C \ R, and hence, a forteriori, for all such z, 

R z (L n )t {anM R Z (L) (n -> 00) 

in the strong operator topology of L 2 (a, b). 15 

Theorem 2 (Stolz and Weidmann (1993, Thm. 4/5)). Let the selfadjoint operators L n 
and L satisfy the assumptions stated above and let Cbe a core ofL such that, eventually, 
CCD(L„). 

(i) lf~L n u^t Lufor all u e C, then L n — > L. 

gsc 

(ii) If L n — > L and if the endpoints of the interval A C R do not belong to the pure 
point spectrum o~pp(L) ofL, the spectral projections to A converge as 

lA(Ln)l( 8 „ A ) ^Ia(L). 



I 4For the cases ^ a < 1 and ^ fi < 1, see the justification of the limit given in Footnote 13. 
I ^We denote by R 2 (L) = (L — z) -1 the resolvent of an operator L. 
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A. 2. Generalized Eigenfunction Expansion of Sturm-Liouville Operators. Let 

L be a formally selfadjoint Sturm-Liouville operator on the interval (a, b), 

Lu = —(pu 1 )' + qu, 

with smooth coefficient functions p > and q. We have the limit point case (LP) 
at the boundary point a if there is some c £ (a, b) and some z 6 C such that 
there exists at least one solution of (L — z)u = Oin {a, b) for which u ^ L 2 (a, c); 
otherwise, we have the limit circle case (LC) at a. According to the Weyl alternative 
(Weidmann 1980, Thm. 8.27) this is then true for all such c and 2 and, moreover, in 
the LP case there exists, for every z £ C \ K, a one-dimensional space of solutions 
u of the equation (L — z)u = for which u £ L 2 (a,c). The same structure and 
notion holds for the boundary point b. 

Theorem 3. Let L be a formally selfadjoint Sturm-Liouville operator on the interval (a, b) 
as defined above. If there is the LP case at a and b, then L is essentially self-adjoint on 
the domain C™(a,b) and, for 2 £ C \ R, the resolvent R Z {L) = (L — z)" 1 of its unique 
selfadjoint extension (which we continue to denote by the letter L) is of the form 

R z (L)cp(x) = * My)<P(y) d v + M*) j x "&(y)<KyWy) ■ 

(A.i) 

Here u a and «j, are the non-vanishing solutions of the equation (L — z)u = 0, uniquely 
determined up to a factor by the conditions u a £ L 1 {a,c) and Uj, £ L 2 (c,b) for some 
c £ (a, b), and W denotes the Wronskian 

W(u a/ u b ) = p(x)(u' a (x)u b (x) - u a (x)u' b (x)), 

which is a constant for x £ (a,b). 

A more general formulation (and a proof) that includes the LC case can be 
found, e.g., in Weidmann (1980, Thm. 8. 22/8. 29). 16 In the following, we write (A.i) 
briefly in the form 

rb 

R z (L)(p(x)= G z {x,y)(p{y)dy. 

J a 

If the imaginary part of the thus defined Green's kernel G z (x, y ) has finite boundary 
values as 2 approaches the real line from above, there is a simple formula for the 
spectral projection associated with L that often applies if the spectrum of L is 
absolutely continuous. 

Theorem 4. (i) Assume that there exits, as e J. 0, the limit 

7r" 1 ImG A+(e (x,i/) K A (x,y), 

locally uniform in x,y £ [a, b ) for each A £ ]R except for some isolated points A for which 
the limit is replaced by 

elmG A+i£ (x,y) 0. 
Then the spectrum is absolutely continuous, o~(L) — cr ac (L), and, for a Borel set A, 

(t A (L)<p,ip) = ( (K x <p,ip)dA (<p,yeC%{a,b)). (A.2) 

J A 



'See Weidmann (1987, pp. 41-42) for a proof that V) is a core if the coefficients are smooth. 



ON THE SCALING LIMITS OF DETERMINANTAL POINT PROCESSES 13 

(ii) Assume further, for some (a' ,b') C (a,b), that 



L L (jj K ^ X '^ d ^) dxd V< co- 



Then 1a(L)1/ 8 / b n is a Hilbert-Schmidt operator on L 2 (a, b) with kernel 



\',V)iy) J K x {x,y)d\, 



If J A Kx(x,y) dX is a continuous function of x,y £ (a', b'), 1^(1)1^;^/) is a trace class 
operator with trace 

trlA(L)l( fl /^/) = J j K\(x,x) dXdx. 

Proof. With £ denoting the spectral resolution of the selfactjoint operator L, we 
observe that, for a given <p £ C™(a,b), the Borel-Stieltjes transform of the positive 
measure }i<p{X) = (E(A)(p,<p) can be simply expressed in terms of the resolvent as 
follows (see Lax 2002, §32.1): 

^ = «,(L)M>. 
, A — z 

If we take z = A + ie and let e \. 0, we obtain by the locally uniform convergence 
of the integral kernel of R z that there exits either the limit 

7r- 1 Im(R A+/e (L)^^) -> (Ktf,<t>) 

or, at isolated points A, 

eJm(R x +ie(L)<p,<l>)^0. 
By a theorem of de la Vallee-Poussin (see Simon 2005, Thm. n.6(ii/iii)), the 
singular part of }iq vanishes, ^^ /S ing = 0; by Plemelj's reconstruction the absolutely 
continuous part satisfies (see Simon 2005, Thm. n.6(iv)) 

d^,ac(A) = (K A <p,<p}dA. 

Since C™{a,b) is dense in L 2 (a,b), approximation shows that E s ing = 0, that is, 
c(L) = <r ac (L) . Since (1a {L)<p, <p) = f A dy.$ (A), we thus get, by the symmetry of the 
bilinear expressions, the representation (A.2), which finishes the proof of (i). The 
Hilbert-Schmidt part of part (ii) follows using the Cauchy-Schwarz inequality and 
Fubini's theorem and yet another density argument; the trace class part follows 
from Gohberg, Goldberg and Krupnik (2000, Thm. IV.8.3) since \r a imJ.^{L)l.r ll im 
is a selfadjoint, positive-semidefinite operator. □ □ 

We apply this theorem to the spectral projections used in the proof of Theorem 1. 
The first two examples could have been dealt with by Fourier techniques (Tao 2012, 
§3.3); applying, however, the same method in all the examples renders the approach 
more systematic. 

Example 1 (Dyson kernel). Consider Lu = —u" on (—00,00). Since u = 1 is a 
solution of Lu — 0, both endpoints are LP; for a given Im 2 > the solutions u a 
(«{,) of (L — z)u = being L 2 at —00 (00) are spanned by 

u a (x)=e- ix ^ z , u b (x)=e ix ^. 
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Thus, Theorem 3 applies: L is essentially selfadjoint on C~(— 00,00), the resolvent 
of its unique selfadjoint extension is represented, for Imz > 0, by the Green's 
kernel 

i f e»'(*-y)V« x>y 

y e -i{x-y)y/z otherwise. 

For A > there is the limit 

-It r ( \ 1/ 1 \ COs((x-y)-v/A) 



2/rVA 

for A < the limit is zero; both limits are locally uniform in x, y 6 ]R. For A = 
there would be divergence, but we obviously have 

e Im Gi e {x, y) — > (e J. 0), 

locally uniform in x,y 6 ]R. Hence, Theorem 4 applies: c(L) = cr ac (L) = [0, 00) 
and (A.2) holds for each Borel set AcR. Given a bounded interval (a, b), we may 
estimate for the specific choice A = (— 00, n 2 ) that 

^ 9 

rb r b 



[fir \K x (x,y)\dX) dxdy 

Ja Ja \ J — 00 / 



b rb I rn 2 
a J a \ JO 



cos((x — y)v A) 



27T\/A 



2 / , , \ 2 



d\\ dxdy <( /'* r -^=J =(&-fl) 2 



« 27iv / A 



T C ° S(( * ^ )VA) rfA= f A sin((x _ !/) yX)dA = sin(7r(x-y)). 
Jo 2vA «A 



Therefore, Theorem 4 yields that \^_ co n 2^{L)\^ a ^ is Hilbert-Schmidt with the 

Dyson kernel 

f 2 K A (*,y) dA = f y)^) AK = M<x-y)) 

J -co v ' Jo 2/rx/A n{x-y) 

restricted to x, y 6 (fl, fr). Here, the last equality is simply obtained from 

i((»- y )VA) ^ _ r 2 rf 

Since the resulting kernel is continuous for x,y E (a,b), Theorem 4 gives that 
l(_co /7T 2) (L)1( S/ ^ is a trace class operator with trace 

trl ( _ 00/7r2) (L)l (fl/b) = fe-«. 

To summarize, we have thus obtained the following lemma. 

Lemma A.i. The operator Lu = —u" is essentially selfadjoint on 00, 00). The 

spectrum of its unique selfadjoint extension is 

a(L) =cr ac {L) = [0,oo). 

Given (a, b) bounded, If^ n i\ {L)1i a ^ is trace class with trace b — a and kernel 

f" 2 cos((x - y)VX) sin(n(x-y)) 
X DyS on(x,y) = / o dX= n{x _ y y . (A.3) 
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Example 2 (Airy kernel). Consider the differential operator Lu = —u" + xu on 
(—00,00). Since the specific solution u(x) = Bi(x) of Lu = is not locally L 2 at 
each of the endpoints, both endpoints are LP. For a given Im z > the solutions 
u a (Ujj) of (L — z)u = being L 2 at —00 (00) are spanned by (see Abramowitz and 
Stegun 1964, Eq. (10.4.59-64)) 

M,i(je) = Ai(x — z) — iBi(x — z), uj,(x) = Ai(x — z). 

Thus, Theorem 3 applies: L is essentially selfadjoint on Cq°( — 00, 00), the resolvent 
of its unique selfadjoint extension is represented, for Imz > 0, by the Green's 
kernel 

Ai(x — z) (Ai(y — z) — z'Bi(y — z)) x > y 

Ai(i/ — z) (Ai(x — z) — z'Bi(x — z)) otherwise. 

For AeR there is thus the limit 

7r _1 ImG A+i0 (^,i/) = K A (x,y) = Ai(x - A)Ai(y - A), 

locally uniform in x, y 6 ]R. Hence, Theorem 4 applies: c(L) = c ac (L) = ]R and 
(A. 2) holds for each Borel set A C R. Given s > —00, we may estimate for the 
specific choice A = (— 00, 0) that 

x 2 \ 

\ *■ \ P CO POO 

\K A (x,y)\dA) dxdy) < / / Ai(x + A) 2 d\dx — t(s) 



G z (x,y) = in 



with 

T (s) = I (2s 2 Ai(s) 2 - 2sAi'(s) 2 - Ai(s)Ai'(s)) . 



3 

Therefore, Theorem 4 yields that l(_oo,o) (^)l(s/») i s Hilbert-Schmidt with the Afn/ 
^° /- 00 Ai(x)Ai'(y) - Ai'(x)Ai(y) 



/U /-oo 
K A (x,y)dA= / Ai(x + A)Ai(y + A) dX 
-00 Jo 



x-y 

restricted to x,y G (s,oo). Here, the last equality is obtained from a Christoffel- 
Darboux type of argument: First, we use the underlying differential equation, 

x Ai{x + A) = Ai" (x + A) - A Ai(x + A), 

and partial integration to obtain 

x / Ai(x + A)Ai(y + A) dX 
Jo 

= / Ai"(x + A)Ai(y + X)dX— / A Ai(x + X)Ai{y + A) dX 
Jo Jo 



r CO r CO 

-Ai'(x)Ai(y) -J Ai'(x + A)Ai'(y + A) dX - X Ai(x + A)Ai(y + A) 



dX. 



Next, we exchange the roles of x and y and substract to get the assertion. Since 
the resulting kernel is continuous, Theorem 4 gives that l^^^^L)!^^ is a trace 
class operator with trace 

trl(_oo,o)( L ) 1 (s / oo) = t(s) ^00 (s-> -00). 
To summarize, we have thus obtained the following lemma. 
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Lemma A.2. The differential operator Lu = —u" + xu is essentially selfadjoint on 
Cg°( — 00,00). The spectrum of its unique selfadjoint extension is 

tr(L) = cr ac (L) = (-00,00). 

Given s > —00, the operator l(_oo,o) (^)l(s,t») *' s ^ race c ' flSS w ^ kernel 

„ , \ r w ^ \ . . , Ai(x)Ai'(y) - Ai'(x)Ai(y) /A , 

K A iry(*,y) = / Ai(x + A)Ai(y + A) dX = — ^ ^— ^ — (A.4) 

Jo x y 

Example 3 (Bessel kernel). Given a > 0, take L« = — 4(xu')' + a 2 x _1 w on (0, 00). 
Since a fundamental system of solutions of Lm = is given by u(x) = x ±a ^ 2 , the 
endpoint x = is LP for a. ^ 1 and LC otherwise; the endpoint x = 00 is LP in 
both cases. Fixing the LP case at x — 0, we restrict ourselves to the case a ^ 1. 

For a given Imz > the solutions u a («/,) of (L — z)« = being L 2 at (00) are 
spanned by (see Abramowitz and Stegun 1964, Eq. (9.1.7-9) and (9.2.5-6)) 



««(*) = Ja(Vxz), U b (x) = Ju(VXz) +iY a (y/xz). 

Thus, Theorem 3 applies: L is essentially selfadjoint on 0^(0, 00), the resolvent of 
its unique selfadjoint extension is represented, for Im z > 0, by the Green's kernel 



. in \Jk(Vxz) (h(y/yz) + iY a (y/yz)) x>y 
G z (x,y) = — < 

4 [}«( Vyz) ( Ja ( Vxz ) + J Y« ( y^z) ) otherwise. 



For A > there is the limit 

1 

7r _1 ImG A+i0 (x,y) = K A (x,y) = -/ a (v / xA)/ a (v^A), 

for A ^ the limit is zero; both limits are locally uniform in x, y 6 ]R. Hence, 
Theorem 4 applies: o~(L) = u ac (L) = [0, 00) and (A. 2) holds for each Borel set 
AcR. Given ^ s < 00, we may estimate for the specific choice A = (—00, 1) that 



1/2 ^ 

Kx{x,y)\dX) dxdy) / / J a (Vx~X) 2 dXdx = T a (s). 

Jo \ ' / / 4 Jo Jo 



Therefore, Theorem 4 yields that l(_oo,i) (^)l(o,s) is Hilbert-Schmidt with the BesseJ 
kernel 

/•1 1 /i — 

/ K A (x,y)d\=- J a (VxX)J a ( v / yX)dX 
J -co 4 Jo 

h(Vx)VyJ' a (Vy) -V^UV^)UVy) 

2(x-y) 

restricted to x, y £ (0, s). Here, the last equality is obtained from a Christoffel- 
Darboux type of argument: First, we use the underlying differential equation, 



X J a (VxA) 



-4~ ( A ^/«(^)) +a 2 A- 1 /«(VxA), 
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and partial integration to obtain 



rfA 



A-^ /«(VxA)) / a (VyA) rfA + ^ jf 1 \- l MVxX)h(y/yX) dX 



Next, we exchange the roles of x and 1/ and substract to get the assertion. Since 
the resulting kernel is continuous, Theorem 4 gives that ^ (L)1(q /S ) is a trace 
class operator with trace 

trl(_oo,i)(L)l( ,s) = T «0) -> 00 (s -> °°)- 
To summarize, we have thus obtained the following lemma. 

Lemma A.3. Given a ^ 1, t/ze differential operator Lu = — 4(x«')' + ot 2 x~^u is essen- 
tially self adjoint on C~(0, 00). The spectrum of its unique self adjoint extension is 

cr{L)=a ac {L) = [0,oo). 

Given ^ s < 00, the operator l(_oo,i) (£)l(o,s) !S ^ race c ' flSS w ^ kernel 

41ei(^y) = J /«(v^A)/«(VyA) rfA 

" 2^ ■ (A ' 5) 

Remark A.i. Lemma A.3 extends to ^ a < 1 if we choose the particular self- 
adjoint realization of L that is defined by the boundary condition (1.6), cf. Hutson, 
Pym and Cloud (2005, Example 10.5.12). 

References 

Abramowitz, M. and Stegun, I. A.: 1964, Handbook of mathematical functions with formulas, graphs, and 

mathematical tables, National Bureau of Standards, Washington. 
Anderson, G. W., Guionnet, A. and Zeitouni, O.: 2010, An Introduction to Random Matrices, Cambridge 

University Press. 

Andrews, G. E., Askey, R. and Roy, R.: 1999, Special Functions, Cambridge University Press, Cambridge. 
Borodin, A. and Olshanski, G.: 2007, Asymptotics of Plancherel-type random partitions, /. Algebra 
313, 40-60. 

Collins, B.: 2005, Product of random projections, Jacobi ensembles and universality problems arising 

from free probability, Probab. Theory Related Fields 133, 315-344. 
Deift, P.: 1999, Orthogonal polynomials and random matrices: a Riemann-Hilbert approach, AMS, Providence. 
Erdelyi, A., Magnus, W., Oberhettinger, F. and Tricomi, F. G: 1953, Higher transcendental functions. Vol. 

II, McGraw-Hill, New York. 
Gohberg, I., Goldberg, S. and Krupnik, N.: 2000, Traces and determinants of linear operators, Birkhauser 

Verlag, Basel. 

Hutson, V., Pym, J. S. and Cloud, M. J.: 2005, Applications of functional analysis and operator theory, second 
edn, Elsevier B. V., Amsterdam. 



i8 



FOLKMAR BORNEMANN 



Johnstone, I. M.: 2001, On the distribution of the largest eigenvalue in principal components analysis, 
Ann. Statist. 29, 295-327. 

Johnstone, I. M.: 2008, Multivariate analysis and Jacobi ensembles: largest eigenvalue, Tracy-Widom 

limits and rates of convergence, Ann. Statist. 36, 2638-2716. 
Kuijlaars, A. B. J.: 2011, Universality, in G. Akemann, J. Baik and P. D. Francesco (eds), Oxford Handbook 

of Random Matrix Theory, Oxford University Press. arXiv : 1103 . 5922. 
Lax, P. D.: 2002, Functional analysis, Wiley-Interscience, New York. 

Simon, B.: 2005, Trace ideals and their applications, second edn, American Mathematical Society, Providence, 
RI. 

Stolz, G. and Weidmann, J.: 1993, Approximation of isolated eigenvalues of ordinary differential 

operators, /. Reine Angeiv. Math. 445, 31-44. 
Tao, T.: 2012, Topics in random matrix theory, American Mathematical Society. 

Wachter, K. W.: 1980, The limiting empirical measure of multiple discriminant ratios, Ann. Statist. 
8, 937-957- 

Weidmann, J.: 1980, Linear operators in Hilbert spaces, Springer- Verlag, New York. 
Weidmann, J.: 1987, Spectral theory of ordinary differential operators, Springer-Verlag, Berlin. 

Zentrum Mathematik - M3, Technische Universitat Munchen, 80290 Munchen, Germany 
E-mail address: bornemannOma . turn . de 



